A conjecture for turbulent flow 
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In this paper, basing on a generalized Newtonian dynamics (GND) approach which has been pro- 
posed elsewhere we present a conjecture for turbulent flow. We firstly utilize the GND to reasonably 
unify the two phenomenological methods recently proposed of the water movement in unsaturated 
soils. Then in the same way a modified Euler equation (MEE) is yielded. Under a zero-order ap- 
proximation, a simple split solution of the MEE can be obtained that shows flow fluids would have 
a velocity field with the power-law scaling feature (power-law fluid) for the case of high Reynolds 
number. 

PACS numbers: 47.27.Ak, 47.27.Jv 



Turbulent flow is one of the most bewildering phe- 
nomena in nature. Ever since a hundred years, people 
have initiated many theories to explain it. These the- 
ories give different phenomenological descriptions, and 
some of them are in good agreement with experimental 
data and provide some very useful clues for our under- 
standing of true essence of turbulent flow. However, up 
to now, we can not employ a unified framework to de- 
scribe turbulent flow yet, not to mention explaining it. 
we even can not give a good definition of turbulence so 
that we can quantitatively determine whether turbulent 
flow has appeared in a fluid. 

In this paper, we would present a conjecture for tur- 
bulent flow with the generalized Newtonian dynamics 
(GND) in order to better describe turbulence. 

The GND describes the fractal world by means of a 
fractional dimensional kinetic velocity (or mass). In the 
anomalous displacement variation model (ADVM), its 
one-dimensional basic dynamics equation in terms of the 
Newtonian kinematics equation form can be written as 
follows: 
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Here m is mass of a particle, Fi is the i-th external force 
acting on the particle, q is a velocity fractal index (vfi) 
and x denotes dx/dt. The right-hand side of the gener- 
alized Newtonian dynamics equation (I) can be under- 
stood as the effective forces acting on the particle in the 
Euclidean space of the fractal environment, and the two 
constants, a and k, which are given by a — m/(qm q ) 
and k = m(l — q) respectively, are the effective force 
coefficients. Here m q is an effective mass of the parti- 
cle in the fractal environment, which has a dimension of 
[kg ■ m 1 ~ q }. The generalized Newtonian dynamics Eq. 
(1) can be reduced to the Newtonian kinematics equa- 
tion mx = £\ Fi when q = I namely the fractal velocity 

±^777- becomes the ordinary velocity 



It is seen that the term k^— in Eq. (1) is an additional 
force which can be eliminated only when the fractal en- 
vironment around the particle disappears, namely q= 1. 

We call the term k^- the fractal force (FF). The FF ap- 
pears in the complex fractal environment, its additional 
property is similar to that of Coriolis force in a rotational 
inertial system, but the origin or the practicable situation 
of it is still unknown for us. 

The cascade process in turbulence shows turbulent flow 
is a type of the fractal phenomenon. Particularly scaling 
law discovered by Kolmogorov p| and the characteristic 
of the j power spectrum of turbulence signal Q can also 
evidently exhibit that. Here we would assume all kinds 
of turbulent flow as well as process from laminar flow 
or convection to turbulence can be described on basis of 
the fractal geometry. The fractal geometry, whose detail 
contents can be reviewed in Ref. J3{, was initialized by 
Mandelbrot. 

Accordingly, we want to introduce the GND into tur- 
bulence. We have a definition of fractal velocity ac- 
cording to the anomalous displacement variation model 
(ADVM): 
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where q is the velocity fractal index (vfi). We assume a 
particle needs a nAt time length to "jump" a (nAx) q dis- 
placement length, then the instantaneous motion state of 
the particle can be expressed by Eq. (2). Here Ax and At 
is a fixed displacement length and time interval respec- 
tively, and n is a natural number representing the n-th 
step of the particle. Because the instantaneous velocity 
of the particle can be written as 
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we let Ax=A\x(t)\, nAx=\x(t)\, and both A|x(t)| and 
At be an arbitrary small increment; allowing for the di- 
rectivity of velocity, Eq. (4) thus can be rewritten with 
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FIG. 1: Schematic representation of the anomalous displace- 
ment variation model (ADVM) and the anomalous time vari- 
ation model (ATVM) of anomalous motion. In the ADVM, 
the particle takes the nAt time length to anomalously move 
(generally jump) the (nAx) q displacement length along the x- 
axis. In the ATVM, the particle needs the (nAt) 9 time length 
to anomalously walk (generally run and wait) the nAx dis- 
placement length along the x-axis. 



x = ±|x| as Eq. (2) where x' = i\x\ q which represents 
the particle's displacement in the complex environment 
(which is somewhat different from real physical displace- 
ment). Similarly, we have another form of definition of 
fractal velocity according to the anomalous time varia- 
tion model (ATVM): 
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In Fig. 1, we give a schematic description of the two 
models. 

But, in fact, Eq. (2) and Eq. (5) are only different 
descriptions of the same physical fact. For example, the 
fact that a particle walks the (nAx) 9 displacement length 
within the nAi time interval can also be stated in another 
word that the particle needs (nAt) p time interval to walk 
nAx displacement length (p ^ q). In addition, \ 1 ~ q x' 
and t are real physical variables in Eq. (2) while x is only 
a scaling variable, where A denotes a characteristic length 
and mA 1_! = m q ; similarly x and r l ~ q t' in Eq. (5) are 
real physical variables while t is only a scaling variable, 
where r is a characteristic time and mr 2q ~ 2 = m'. Note 
that x— space is changed to x 9 — space when a particle 
is moving with time, namely ordinary scaling matching 
structure of time and space is modified, but any force 
(except those related to velocity) is invariant from the 
ordinary environment to the fractal space. 

With the help of the GND approach, we can easily 
unify the two phenomenological methods recently pro- 
posed of the water percolation in unsaturated soils. In 
the unsaturated soil- water transport, the water flux for 
the horizontal one-dimensional column case follows the 
Buckingham-Darcy law: 



Q = -k(0) 



(6) 




where Q denotes the volume flow rate per unit area, k 
is unsaturated hydraulic conductivity, ij) is the hydraulic 
pressure head and 9 is the soil water content; and the 
customary Richards' equation 
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can be arrived at by combining Eq. (6) and the mass 
conservation equation ff + |§ = with D{9) = k{9)^. 
Here D(9) is the soil water diffusivity. 

It is necessary to emphasize that movement of liquid 
water in soils is not a diffusion phenomenon Con- 
versely, it is a kind of macroscopic flow motion which is 
based on Navier-Stokes law. Phillip derived the conclu- 
sion that Buckingham-Darcy law follows from the Navier- 
Stokes equation with (u-'V)u ~ 0: «fcV 2 u = V</> where Vk 
is dynamical viscosity and <fi is total potential Here 
no concept of probability appears. 

However, deviations from Eq. (7) cause further as- 
sumptions that the diffusivity has a dependence on dis- 
tance given by |x| - ^.D(0) or the diffusivity is a time- 
dependent quantity given by t 2m ~ 1 D(8) 0. Here m is a 
positive constant. We now want to unify the above two 
assumptions and do not need to introduce other models 
for microscopic anomalous transport. 

In the GND framework, the reduced Navier-Stokes 
equation is rewritten as v' k V 2 u s = where u s is fractal 
velocity. Then the Richards' equation is generalized as 



and 
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respectively. Here s is the vfi and D s {9) is a fractal soil 
water diffusivity. 

Eq. (8) and Eq. (9) are similar to the generalized 
Richards' equations proposed in [9] and [10], respectively, 
which have been shown to be valid for some experiment 
data. But we can unify the two kinds of assumption 
in one framework, and it seems more reasonable be- 
cause previous frameworks of anomalous diffusion are not 
appropriated for describing essentially macroscopic flow 
movements. 
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For the turbulence problem, we can obtain a modified 
Euler equation (MEE) from the GNKE Eq. (1) according 
to the principle of the GND: 

dv o. ol -* 
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FIG. 2: Velocity field v x (x, y) given by a solution of the MEE 
Eq. (19) for different values of q. The solution reflects motion 
state of a flow in the fractal environment. Here a and b are 
both taken 1. 



p is density of a fluid; —adW q p is an effective pressure 
acting on certain a volume element dV from surrounding 
fluid bodies and F q is considered as an effective complex 
force field acting on a volume element dV of the fluid 
where dFi (i = x, y, z) is a real physical force (other than 
pressure) acting on a volume element dV . We notice that 
df has the same dimension as another nonlinear term 
(v-V)v. 

We call V g the fractional gradient operator, and also 
we notice that pf3df is exactly the fractal force on a unit 
volume of fluid. In addition, the coefficient a with the 
dimension [m 9_1 ] can be seen as A 9-1 . Here A is a char- 
acteristic length of the fluid. We see the MEE (10) will 
will recover the ordinary Euler equation: 

+ (ff.y )v = — Vp+-— . 14 
at p p dv 

Here dF is force (other than pressure) acting on the vol- 
ume element dV . 

The MEE reflects the situation of a fluid with fractal 
feature, since it is derived directly from a fractal velocity 
which describes the fractal state of a system. We think 
the ordinary Euler equation can not play a role in de- 
termining the fractal state of the fluid, for the nonlinear 
term (v ■ \7)v is only a kinematic factor rather than a 
dynamical factor like the term (1 — q)df in the MEE. In 
other words, when the fluid becomes in a fractal state, 
the fluid would be acted on by a special fractal force 
which appears only in the fractal situation, and the frac- 
tal force can also be seen as an external energy input; 
at the same time, ordinary forces will also be deformed, 
but these changes can not follow naturally from the or- 
dinary Euler equation. Of course, the process from a 
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laminar flow or a convection current to turbulent flow is 
controlled by real physical parameters such as velocity of 
the fluid or temperature difference, these parameters of- 
ten can be combined into a dimensionlcss number such as 
the Reynolds number Re. Particularly there is a critical 
Reynolds number Rec. When the Reynolds number of a 
fluid is above Rec, turbulence happens. That can make 
us guess that the vfi q is a function of these dimensionless 
numbers and critical dimensionless numbers such as Re 
and Rec such that the fact that turbulence emerges when 
Re > Rec is equal to the situation of the MEE for q > f 
(or q < 1). 

The viscous term of a fluid can be directly added in F q 
in Eq. (10), but here we would take a look at solutions 
of the MEE in a range of high Reynolds number. 

We study the MEE for the case of the two-dimensional 
incompressible and constant flow, and we consider F q — 
0; then Eq. (10) becomes 



v,\yr 1 ^+v y \y\^ 9 -^ = -^+p\ y r^ 



We make such an approximation that 
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where a and b, which represent the mean velocity of the 
fluid v x and v y , respectively, are constant. The approxi- 
mation is reasonable for a fluid with small velocity gra- 
dient. 

Eliminating the pressure intensity p and utilizing the 
continuity condition, 



dv x dv y 
dx dy 



we have 



■'0 + H-0 = 
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where ip = bv x 



av v 



We let vi = 



2 _ (jj+avy) 



— ~ ip 2 as well as v 2 
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ip 2 , which is 



proper for a fluid whose velocity is not too large. Then we 
obtain a second-order partial differential equation con- 
cerning ip: 
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dx 
(19) 

We neglect the two second partial derivative terms of 
ip, because both of them can be seen as small quanti- 
ties. Eq. (19) thus is simplified as a first-order partial 
differential equation: 
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whose general solution can be easily arrived at: 

iP = *(«) = ^{a 2 ^- 1 + b 2 ^- 1 ) . (21) 

Here $ is any continuous and differentiable function of 
u. 

For simplicity, we choose ip as follows: 

bv x - av y = (b- aX^laf- 1 + b 2 ^ 1 ) , (22) 
such that 



v« = a 2 ^- 1 



b 2 \y\ q ~ 



(23) 



Evidently, the velocity field v x or v y is a parabolic dish 
for different values of q. 

Therefore, under such an approximation, we obtain a 
simple solution of the MEE for a two-dimensional incom- 
pressible and constant flow in the high Reynolds number 
limit: 



vJ\Xi 



(24) 



where v^j denotes Uj or Vj. 

We graph velocity field v x given by Eq. (23) for dif- 
ferent values of q in Fig. 2. Here a and b are both taken 
1. It is clearly seen that manifold is broken when q < 1 
and the fluid velocity v x becomes infinity at the position 
x = or y = 0. We can imagine that the fluid body 
could be broken up into four closed packets in each di- 
mension if proper boundary conditions about (x,y,Vij) 
are given. Thus we can call the solution Eq. (23) for the 
case of q < 1 a split solution for the MEE. The situation 
for 1 < q < 2 is somewhat different where four convexes 
of the fluid are not broken. When q > 2, velocity field of 
the fluid becomes an concave parabolic dish. For q = 2 
the velocity field comprises four planes and for q = 1 a 
laminar flow case is recovered. 

The solution reflects motion state of a flow in the frac- 
tal environment which is embodied in the MEE by both 
a fractal force which can be understood as an external 
nonlinear energy input and a fractional gradient opera- 
tor V g which is similar to the Riesz fractional derivative 
operator ||| ^| . It is important that the fractional gradi- 
ent operator will simultaneously be reduced to the ordi- 
nary gradient operator when the fractal force disappears. 
However, whether the MEE can describe turbulent flow 
or not is still an open problem. 

Summarizing, the fact that the GND framework can 
reasonably unify the two phenomenological methods re- 
cently proposed of anomalous transport of water in un- 
saturated soils and other satisfactory results of the GND 
elsewhere make us conjecture that the GND constitutes 
a dynamical basis of the problem of turbulence, thus a 
modified Euler equation (MEE) is yielded. Under a zero- 
order approximation, a simple split solution of the MEE 
can be obtained, and we see that turbulent flow would 
have a velocity field with the power-law scaling feature 
for the case of high Reynolds number. 
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